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Abstract 

We consider a quantum system consisting of a regular chain of elementary subsystems with nearest neighbor interactions 
and assume that the total system is in a canonical state with temperature T. We analyze under what condition the state 
factors into a product of canonical density matrices with respect to groups of n subsystems each, and when these groups 
have the same temperature T. While in classical mechanics the validity of this procedure only depends on the size of the 
groups n, in quantum mechanics the minimum group size n m - ln also depends on the temperature T ! As examples, we apply 
our analysis to a harmonic chain and an Ising spin chain. We discuss various features that show up due to the characteristics 
of the models considered. For the harmonic chain, which successfully describes thermal properties of insulating solids, our 
approach gives a first quantitative estimate of the minimal length scale on which temperature can exist: This length scale 
is found to be constant for temperatures above the Debye temperature and proportional to T" 3 below. 
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1. Introduction 

The microscopic limit of the applicability of Ther- 
modynamics is still poorly understood [1,2]. Down to 
which length scales can its standard concepts mean- 
ingfully be defined and employed? However, despite its 
fundamental relevance this topic has only received lit- 
tle attention for a long time [9,10]. 

The situation drastically changed in recent years, 
since the issue became highly relevant for experiments. 
Significant progress in the synthesis and processing 
of materials with structures on nanometer length 
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scales has created a demand for better understanding 
of thermal properties of nanoscale devices, individ- 
ual nanostructures and nanostructured materials [3]. 
Experimental techniques have improved to such an 
extent that the measurement of thermodynamic quan- 
tities like temperature with a spatial resolution on the 
nanometer scale seems within reach [4,5,6,7]. 

To provide a basis for the interpretation of present 
day and future experiments in nanoscale physics and 
technology and to obtain a better understanding of the 
limits of thermodynamics, it is thus indispensable to 
clarify the applicability of thermodynamical concepts 
on small length scales starting from the most funda- 
mental theory at hand, i. e. quantum mechanics. 
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Here we focus on the question whether tempera- 
ture can be meaningfully defined on those small length 
scales. 

The existence of thermodynamical quantities [11], i. 
e. the existence of the thermodynamic limit strongly 
depends on the correlations between the considered 
parts of a system. Hence, local temperature is only ex- 
pected to exist in systems with short range interac- 
tions, where the interaction between a region in space 
and its surrounding only takes place at the surface of 
that region. Being only a surface effect, correlations 
between the region and the surrounding thus become 
less important as the size of the region increases, since 
the energy contained in the region grows with its vol- 
ume (diameter ) while correlations only grow with the 
surface (diameter 2 ) [8]. 

This scaling of interactions between parts of a sys- 
tem [12] thus sets a minimal length scale on which cor- 
relations are still small enough to permit the definition 
of local temperatures. In this paper we study this con- 
nection quantitatively. 

We adopt here the convention that a local tempera- 
ture exists if the considered part of the system is in a 
canonical state, where the distribution is an exponen- 
tially decaying function of energy characterized by one 
single parameter. This implies that there is a one-to- 
one mapping between temperature and the expectation 
values of observables, by which temperature is usually 
measured. Temperature measurements based on differ- 
ent observables will thus yield the same result, contrary 
to distributions with several parameters [24] . In large 
systems composed of very many subsystems, the den- 
sity of states is a strongly growing function of energy 
[14]. If the distribution were not exponentially decay- 
ing, the product of the density of states times the dis- 
tribution would not have a pronounced peak and thus 
physical quantities like energy would not have "sharp" 
values. 

Based on the above arguments and noting that 
a quantum description becomes imperative at 
nanoscopic scales, the following approach appears to 
be reasonable: Consider a large homogeneous quan- 
tum system, brought into a thermal state via inter- 
action with its environment, divide this system into 
subgroups and analyze for what subgroup-size the 
concept of temperature is still applicable. 

Harmonic lattice models are a standard tool for the 
description of thermal properties of solids. We there- 



fore apply our theory to a harmonic chain model to 
get estimates that are expected to be relevant for real 
materials and might be tested by experiments. 

Recently, spin chains have been subject of extensive 
studies in condensed matter physics and quantum in- 
formation theory. Thus correlations and possible local 
temperatures in spin chains are of interest, both from a 
theoretical and experimental point of view [13,16,15]. 
We study spin chains with respect to our present pur- 
pose and compare their characteristics with the har- 
monic chain. 



2. General Theory 

We consider a homogeneous (i.e. translation invari- 
ant) chain of elementary quantum subsystems with 
nearest neighbor interactions. The Hamiltonian of our 
system is thus of the form 

H = Y / Hi + I i} i+i, (1) 

i 

where the index i labels the elementary subsystems. 
Hi is the Hamiltonian of subsystem i and the 
interaction between subsystem i and i + 1. We assume 
periodic boundary conditions. 

We now form Ng groups of n subsystems each (index 
i — » (n — l)n + j; fi = 1, . . . , Ng; j = 1, . . . , n) and 
split this Hamiltonian into two parts, 

H = H + I, (2) 

where Ho is the sum of the Hamiltonians of the isolated 
groups, 

N G 

H = ^2 (Hp ~ witn 

n 

H-V- = yi Hn(n-l)+j + In(n-l)+j, n( M -l)+j + l (3) 
J = l 

and I contains the interaction terms of each group with 
its neighbor group, 

N G 

I = ^ (4) 

f»=l 

The eigenstates of the Hamiltonian Ho are products of 
group eigenstates, 

N G 

H | a) = E a |a) with |o> = |o M >, (5) 
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where (7Y M - J M „, M „+i) |o M ) = £ M ja M ). £ M istheenergy 
of one subgroup only and E a = YLu=i 



2.1. Thermal State in the Product Basis 



We assume that the total system is in a thermal state 



(6) 



where Z is the partition sum and [3 = (fcsT) -1 the 
inverse temperature with Boltzmann's constant fcg and 
temperature T. 

We are interested in the matrix representation of the 
state (6) in the product basis, the eigenbasis of Ho. 
The diagonal elements are the expectation values of 
the density operator (6) in the states \a), 



r-E x 
(a\p\a) = I 
Je q 



Wa(E) 



-0E 



dE, 



(7) 



where Eo is the energy of the ground state and E\ the 
upper limit of the spectrum [19]. w a (E) is defined as 
follows: The state \a) is not an eigenstate of the to- 
tal Hamiltonian H. Thus, if H would be measured in 
the state \a), eigenvalues of H would be obtained with 
certain probabilities: w a (E) is the density of this prob- 
ability distribution. Since the Hamiltonian H is the 
sum of Hamiltonians of the groups, the situation has 
some analogies to a sum of random variables. This in- 
dicates that there might exist a central limit theorem 
for the present quantum system, provided the number 
of groups becomes very large [18] . Since the state \a) 
is not translation invariant and since H also contains 
the group interactions, the central limit theorem has 
to be of a Lyapunov (or Lindeberg) type for mixing 
sequences [17]. One can indeed show that such a quan- 
tum central limit theorem exists for the present model 
[20,23] and that w a (E) thus converges to a Gaussian 
normal distribution in the limit of infinite number of 
groups N G , 



lim w a {E) 

Nq — *oo 



27rA a 



exp 



(E — Eg — Eg) 

2A2 



where the quantities e a and A a are defined by 

e a = (a\H\a) - E a 
A 2 a = (a\H 2 \a}-(a\H\a} 2 . 



(8) 

(9) 
(10) 



e a is the difference between the energy expectation 
value of the distribution w a (E) and the energy E a , 



while A 2 is the variance of the energy E for the distri- 
bution w a (E). Note that E a has a classical counterpart 
while A 2 is purely quantum mechanical. It appears be- 
cause the commutator [H, H ] is nonzero, and the dis- 
tribution w a (E) therefore has nonzero width. The two 
quantities s a and A^ can also be expressed in terms of 
the interaction only (see eq. (2)), 



E a = (a|/|o) 

A 2 a = (a\I 2 \a)-(a\I\a) 2 



(11) 
(12) 



meaning that £ a is the expectation value and A 2 the 
squared width of the interactions in the state \a). 

The rigorous proof of equation (8) is given in [20] and 
based on the following two assumptions: The energy of 
each group H M as defined in equation (3) is bounded, 
i. e. 

(x\nM < c (13) 

for all normalized states |x) and some constant C, and 
(a\H 2 \a) - (a\H\a) 2 > N G C (14) 

for some constant C' > 0. 

In scenarios where the energy spectrum of each ele- 
mentary subsystem has an upper limit, such as spins, 
condition (13) is met a priori. 

For subsystems with an infinite energy spectrum, 
such as harmonic oscillators, we restrict our analysis 
to states where the energy of every group, including 
the interactions with its neighbors, is bounded. Thus, 
our considerations do not apply to product states \a), 
for which all the energy was located in only one group 
or only a small number of groups. Since Nq S> 1, the 
number of such states is vanishingly small compared 
to the number of all product states. 

If conditions (13) and (14) are met, equation (7) can 
be computed for iV G > 1 [21]: 



(o|p|o) 



erfc 



1 / — B 2 A 2 



E Q -E a + !3A 2 g 
V2A a 



erfc 



E!-E a + PA 2 a 



V2A a 



(15) 



where E a = (a\H\a) = E a + £ a and erfc(x) is the con- 
jugate Gaussian error function [25]. The second error 
function in (15) only appears if the energy is bounded 
and the integration extends from the energy of the 
ground state Eo to the upper limit of the spectrum Ei . 
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Note that E a is a sum of Ng terms and that A a 
fulfills equation (14). The arguments of the conju- 
gate error functions thus grow proportional to \JNg 
or stronger. If these arguments divided by \/Ng are 
finite (different from zero), the asymptotic expan- 
sion of the error function [25] may thus be used for 
Ng 3> 1: erfc(x) « exp (— x 2 ) / (^/tyx) for x — > oo 
and erfc(x) ~ 2 + exp (— x 2 ) / (^/tyx) for x — > — oo. 
Inserting this approximation into equation (15) and 
using Eo < E a < Ei shows that the second conjugate 
error function, which contains the upper limit of the 
energy spectrum, can always be neglected compared 
to the first, which contains the ground state energy. 

The off diagonal elements {a\p\b) vanish for 
\E a — Et\ > A a + A;, because the overlap of the two 
distributions of conditional probabilities becomes neg- 
ligible. For \E a — Eb\ < A a + Ab, the transformation 
involves an integral over frequencies and thus these 
terms are significantly smaller than the entries on the 
diagonal. 

2.2. Conditions for Local Thermal States 

We now test under what conditions the density ma- 
trix p may be approximated by a product of canonical 
density matrices with temperature /3i oc for each sub- 
group fj, = 1,2,..., Ng- Since the trace of a matrix is 
invariant under basis transformations, it is sufficient 
to verify the correct energy dependence of the prod- 
uct density matrix. If we assume periodic boundary 
conditions, all reduced density matrices are equal and 
their product is of the form (a\p\a) oc cxp(— P\ oc E a ). 
We thus have to verify whether the logarithm of rhs of 
equation (15) is a linear function of the energy E a , 

ln((a\p\a))z>-0i oc E a + c, (16) 

where 0\ oc and c are constants. 

Applying the asymptotic expansion of the error func- 
tion in equation (15), one realizes that equation (16) 
can only be true for 



Ea+e a 



'Ng A a 



and 



+ ^A 2 
+ 2 a 



^NaAa 



C\E a + C 2 , 



(17) 



(18) 



where ci and C2 are arbitrary real constants. Note that 
e a and A 2 need not be functions of E a and therefore 
in general cannot be expanded in a Taylor series. 



To ensure that the density matrix of each subgroup 
p is approximately canonical, one needs to satisfy (18) 
for each subgroup p separately; 

£ M -1 +S U , 2 A 2x 72 



Cl Ep, + C2 

(19) 

wheree M = (a|/ Mn , M „ + i|a), A 2 = {a\H 2 u \a) - {a\H^\a} 2 
and A 2 = <a|«„-i«„ + HMv-i\a)- 

-2££i-iM«,-i|a><a|«,,|a). 

Temperature becomes intensive, if the constant ci 
vanishes, 

|ci|<l /3ioc=/3. (20) 

If this was not the case, temperature would not be 
intensive, although it might exist locally. 

It is sufficient to satisfy conditions (17) and (19) for 
an adequate energy range _B m in < E^ < iJmax only. 
The density of states of large many body systems is 
typically a rapidly growing function of energy [26,14]. If 
the total system is in a thermal state, occupation prob- 
abilities decay exponentially with energy. The product 
of these two functions is thus sharply peaked at the 
expectation value of the energy E of the total system 
E + Eo =Tr(Hp), with Eo being the ground state en- 
ergy. The energy range needs to be centered around 
this peak and large enough. On the other hand it must 
not be larger than the range of values E u can take on. 
Therefore a pertinent and "safe" choice for i? m in and 



Smin = max([^] min ,iA + ^ 

L '"<>« - ' Jml ( [E ^' a JtG~ + lt 



(21) 



where a ^> 1 and E will in general depend on the global 
temperature. In equation (21), [-E M ] min and [£ M ] max de- 
note the minimal and maximal values Ea can take on. 

For a model obeying equations (13) and (14), the two 
conditions (17) and (19), which constitute the general 
result of this article, must both be satisfied. These fun- 
damental criteria will now be applied to some concrete 
examples. 



3. Results for Special Models 

We are now going to examine equations (17) and 
(19) for some special models. As discussed in the intro- 
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duction, those terms in (17) and (19), which contain 
the interactions between neighboring groups (e M and 
A M ) are independent of the group size, while the "back- 
ground" term E^ grows approximately linear with the 
group size n. Hence, the interaction terms become less 
important compared to with increasing group size. 
Therefore, taking into account the energy range (21), 
conditions (17) and (19) can be used to determine a 
minimal group size n m i n . Technical details of this cal- 
culation can be found in [21] and [22]. 

For each model, the results depend, besides the 
global temperature /3, on two "accuracy" parameters, 
a and 8, which quantify the tolerated deviations from 
a canonical distribution. The first parameter, a, has 
already been introduced in equation (21). It is a mea- 
sure for the energy range, where one demands (17) 
and (19) to be fulfilled. The second parameter, 8, is 
the ratio between the terms in the lhs of condition 
(19), which are linear in E^, and those, which are of 
higher order. Thus 8 — 0.01 means that all terms in 
the lhs of (19) which are not linear in _B M are at least 
two orders of magnitude smaller than the linear ones. 



3.1. Harmonic Chain 

As a representative for the class of systems with 
an infinite energy spectrum, we consider a harmonic 
chain of Ng ■ n particles of mass m and spring constant 
y^rnujo- In this case, the Hamiltonian reads 

Hi = — Pi + — w 9, (22) 
h,i+i = -muj q l q i+1 , (23) 

where p» is the momentum of the particle at site i and 
qi the displacement from its equilibrium position i ■ ao 
with ao being the distance between neighboring parti- 
cles at equilibrium. We divide the chain into Ng groups 
of n particles each and thus get a partition of the type 
considered above. Describing the groups in Debye ap- 
proximation (see [21] for details) we calculate the min- 
imal group size n m in. 

Here, we only present the results for the present 
model, which are shown in figure 1. 

For high (low) temperatures be estimated 

by 
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Fig. 1. Log-log-plot of n m j n from eq. (17) (dashed line) and 
n min from eq. (19) (solid line) for a = 10 and <5 = 0.01 
as a function of T/0 for a harmonic chain. is the De- 
bye temperature, a characteristic constant of the respective 
material, which can be found tabulated [27]. <5 and a are 
explained at the beginning of section 3. Local temperature 
exists in the shaded area. 

{2a/S forT>9 
(24) 
(3a/27r 2 ) (9/T) 3 for T < 9 

where O is the Debye temperature, a characteristic con- 
stant of the respective material, which can be found 
tabulated [27]. 

As an extra result, we get, that whenever local tem- 
perature exists, it is equal to the global one, i. e. tem- 
perature is intensive. 

3.2. Ising Spin Chain in a Transverse Field 

In this section we consider an Ising spin chain in a 
transverse field. For this model the Hamiltonian reads, 

Hi = -Bat 

h,i+i = (°f ® fff+i ~ o\ ® crf+i) , (25) 

where o~f,a^ and a\ are the Pauli matrices. B is the 
magnetic field and J a coupling parameter. We will 
always assume B > 0. 

Running through a similar calculation as for the har- 
monic chain, the details of which can be found in [22], 
we get the result shown in figure 2. 

As for the harmonic chain, temperature is intensive 
whenever it exists locally, i.e. then local temperatures 
are equal to global ones. 
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Fig. 2. Log-log-plot of n m i n for J = 0.1 X B from eq. (17) 
(dashed line), and from eq. (19) (dash - dotted line) and 
™min f° r J = 10 X B from eq. (17) (solid line) and from eq. 
(19) (gray line) as a function of T/B. a = 10 and <5 = 0.01. 
<5 and a are explained at the beginning of section 3. 

4. Summary and Conclusions 

We have considered a linear chain of particles in- 
teracting with their nearest neighbors. We have par- 
titioned the chain into identical groups of n adjoining 
particles each. Taking the number of such groups to 
be very large and assuming the total system to be in a 
thermal state with temperature T we have found condi- 
tions (equations (17) and (19)), which ensure that each 
group is approximately in a thermal state. Further- 
more, we have determined when the isolated groups 
have the same temperature T, that is when tempera- 
ture is intensive. 

The result shows that, in the quantum regime, these 
conditions depend on the temperature T, contrary to 
the classical case. The characteristics of the tempera- 
ture dependence are determined by the width A a of 
the distribution of the total energy eigenvalues in a 
product state, which in turn is related to the fact that 
the commutator [H, Ho] is nonzero. In particular, the 
ground state of the total system does not factorize with 
respect to any partition [16,29,28,30]. 

We have then applied the general method to a har- 
monic chain and an Ising spin chain in a magnetic field. 
When applied to specific models, conditions (17) and 
(19) determine a minimal group size n m i n for which a lo- 
cal canonical state and hence local temperature can ex- 
ist. We have given order of magnitude estimates of n m in 
for these models. Grains of a smaller size are no more in 
a thermal state and temperature measurements with a 
higher resolution should therefore no longer be inter- 



preted in a standard way. 

For the spin chain and the harmonic chain, the tem- 
perature dependencies of ^min for low temperatures co- 
incide, n m i n oc T~ 3 , because both couplings have the 
same structure and the upper limit of the spectrum of 
the spin chain becomes irrelevant at low temperatures. 

On the other hand, the high temperature asymp- 
totics of n m i n differs significantly: 

For spins at very high global temperatures, the to- 
tal density matrix is almost completely mixed, i. e. 
proportional to the identity matrix, and thus does not 
change under basis transformations. There are thus 
global temperatures which are high enough, so that lo- 
cal temperatures exist even for single spins. 

For the harmonic chain, this feature does not appear, 
since the size of the relevant energy range increases 
indefinitely with growing global temperature, leading 
to the constant minimal length scale in the high energy 
range. 

Since harmonic lattice models in Debye approxima- 
tion have proven to be successful in modeling thermal 
properties of insulators (e.g. heat capacity) [27], our 
calculation for the harmonic chain provides a first es- 
timate of the minimal length scale on which intensive 
temperatures exist in insulating solids, 

'min = flminflO- (26) 

Let us give some numerical estimates: Choosing the 
"accuracy parameters" to be a = 10 and 8 = 0.01, we 
get for hot iron (T > 6 « 470 K, a « 2.5 A) Z min « 
50 fim, while for carbon (O w 2230 K, ao ~ 1.5 A) at 
room temperature (270 K) Z m i n ~ 10 fim. The coarse- 
graining will experimentally be most relevant at very 
low temperatures, where Z m in may even become macro- 
scopic. A pertinent example is silicon (6 w 645 K, ao £s 
2.4 A), which has Z m i n « 10cm at T w 1 K (again with 
ol = 10 and 5 = 0.01). 

Of course the validity of the harmonic lattice model 
will eventually break down at finite, high temperatures 
and our estimates will thus no longer apply there. 

Measurable consequences of the local breakdown of 
the concept of temperature are interesting questions 
which arise in the context of the present discussion. 

In the scenarios of global equilibrium, which we con- 
sider here, a temperature measurement with a micro- 
scopic thermometer, locally in thermal contact with 
the large chain, would not reveal the non existence of 
local temperature. One can model such a measurement 
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with a small system, representing the thermometer, 
coupled to a heat bath, representing the chain. It is 
a known result of such system bath models [31], that 
the system always relaxes to a thermal state with the 
global temperature of the bath, no matter how local 
the coupling might be. 

This, however, does not mean that the existence or 
non existence of local temperatures had no physical rel- 
evance: There are indeed physical properties, which are 
determined by the local states rather than the global 
ones. Whether these properties are of thermal char- 
acter depends on the existence of local temperatures. 
A detailed discussion of such properties will be given 
elsewhere [24]. 

The length scales, calculated in this paper, should 
also constrain the way one can meaningfully define 
temperature profiles in non-equilibrium scenarios [32] . 
Here, temperature measurements with a microscopic 
thermometer, which is locally in thermal contact with 
the sample, might indeed be suitable to measure the 
local temperature. An explicit study of this possibility 
should be subject of future research. 

We thank M. Michel, M. Henrich, H. Schmidt, M. 
Stollsteimer, F. Tonner and C. Kostoglou for fruitful 
discussions. 
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